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Abstract 

We study the controllability of the Bloch equation, for an ensemble of non inter- 
acting half-spins, in a static magnetic field, with dispersion in the Larmor frequency. 
This system may be seen as a prototype for infinite dimensional bilinear systems with 
continuous spectrum, whose controllability is not well understood. We provide several 
mathematical answers, with discrimination between approximate and exact controlla- 
bility, and between finite time or infinite time controllability: this system is not exactly 
controllable in finite time T with bounded controls in L 2 (0, T), but it is approximately 
controllable in L°° in finite time with unbounded controls in Lf£ c ([0, +00)). Moreover, 
we propose explicit controls realizing the asymptotic exact controllability to a uniform 
state of spin +1/2 or —1/2. 

Key words, bilinear control systems, Bloch equation, continuous spectrum, controlla- 
bility of infinite dimensional systems, ensemble controllability, quantum systems. 

1 Introduction 

1.1 Studied system, bibliography 

Most controllability results available for infinite dimensional systems are related to systems 
with discrete spectra. As far as we know, very few controllability studies consider systems 
admitting a continuous part in their spectra. In [12] an approximate controllability result 
is given for a system with mixed discrete/continuous spectrum: the Schrodinger partial 
differential equation of a quantum particle in an N-dimensional decaying potential is shown 
to be approximately controllable (in infinite time) to the ground bounded state when the 
initial state is a linear superposition of bounded states. 

In O [lOl QT] a controllability notion, called ensemble controllability, is introduced and 
discussed for quantum systems described by a family of ordinary differential equations (Bloch 
equations) depending continuously on a finite number of scalar parameters and with a finite 
number of control inputs. Ensemble controllability means that it is possible to find open- 
loop controls that compensate for the dispersion in these scalar parameters: the goal is to 
simultaneously steer a continuum of systems between states of interest with the same control 
input. The articles [lOl |TT] highlight, for three common dispersions in NMR spectroscopy, 
the role of Lie algebras and non-commutativity in the design of a compensating control 
sequence and consequently in the characterization of ensemble controllability. 
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Such continuous family of ordinary differential systems sharing the same control inputs 
can be seen as the prototype of infinite dimensional systems with purely continuous spectra. 
The goal of this paper is to show that the very interesting controllability analysis of [9j[l0j[TT] 
can be completed by functional analysis methods developed for infinite dimensional systems 
governed by partial differential equations (see, e.g., [7] for samples of these methods). 

We focus here on one of the three dispersions cases treated in [SHUHE]- We consider an 






ensemble of non interacting half-spins in a static field [ o ] in M. 3 , subject to a transverse 

V Bo J 

radio frequency field -«(t) in R 3 (the control input). The ensemble of half-spins is 
V o J 

described by the magnetization vector MeK 3 depending on time t but also on the Larmor 
frequency to — —jB (7 is the gyromagnetic ratio). It obeys to the Bloch equation: 

dM / -« v(t) \ 

— (t,w)= u -u(t) \M(t,u), (t,u/)e[0,+oo)x(a;„u>*), (1) 

m V -»(*) «(*) / 

where —00 ^ lu« < uj* ^ +00 are given . With the notations 

/ \ / 1 \ 

Q x := -1 , Cl y := ) , ft* := I 1 I , (2) 



the system ([T]) can be written 

^-(t,Lj) = (wft z +u(t)n x + v{t)n y ) M(t,cu), (t,Lj) G [0, +00) x (u>*,u)*). (3) 

It is a bilinear control system in which, at time t, 

• the state is (M(t,oj))^ e ^ r u ,y, for each cu, M(t,co) € § 2 , the unit sphere of M 3 , 

• the two control inputs u(t) and v(t) are real. 

In the sequel, we denote by ej,, the R 3 -vector of coordinates (Ski)ie{i,2,3}- Thus, we study the 
simultaneous controllability of a continuum of ordinary differential equations, with respect 
to a parameter uo that belongs to an interval Notice that, when v = u = 0, the 

spectrum of this system is made by the union of the two segments, i(w*,w*) and — 
belonging to the imaginary axis. 

The pioneer articles |{J [101 QT| provide convincing arguments indicating why the system 
([3]) is ensemble controllable (i.e. approximately controllable in L 2 ((uj«, u>*), S 2 )) with un- 
bounded and also bounded controls, when and u>* are finite. Here, we provide several 
mathematical results that complete these ensemble controllability results with discrimi- 
nations between approximate or exact controllability and between finite or infinite time 
(asymptotically) controllability. 

1.2 Controllability issues 

Let us recall a famous non controllability result for infinite dimensional bilinear systems due 
to Ball, Marsden and Slemrod [Q. This result concerns general systems of the form 

^-=Aw+p(t)Bw (4) 
where the state is w and the control is p : [0, T] — > ML 
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Theorem 1 Let X be a Banach space with dim(X) = +00, A generate a C° -semigroup of 
bounded operators on X and B : X — > X be a bounded operator. For wq € X, w(t;p,wo) 



is contained in a countable union of compact subsets of X and, in particular, it has an empty 
interior in X. Thus ^ is not controllable in X with controls in U r >iL[ oc ([0, +00)). 

We cannot apply directly here this result since the spaces X = L 2 ((uj*,u}*),§ 2 ) or 
C°([ui*, uj*}, S 2 ) where the Cauchy problem is well-defined are not vector spaces. In or- 
der to get an interesting result for the Bloch equation, one needs extensions of the above 
result to Banach manifolds. (This has been done in [14] when the manifold is the unit sphere 
of a Hilbert space.) For J5J|, the situation is similar to the one described in Theorem [TJ In 
TheoremlU we show that for any analytic initial condition M (uj), the reachable set in finite 
time T > from M with controls in L 2 (0,T) only contains analytic functions of u>, Thus, 
the reachable set (from an analytic initial dara) has an empty interior in L 2 ((w*, to*), § 2 ), 
which is a natural space for the Cauchy problem. 

However, for (J5J), the obstruction to exact controllability given by Theorem [5] has much 
stronger consequence than the obstruction described by Theorem Q] which is, in fact, a rather 
weak non controllability result. Indeed, it does not prevent the reachable set from being 
dense in X (approximate controllability in X). For example, this is the case for the ID 
beam equation 



in which the state is (u,Ut) and the control is p. Theorem Q] ensures that this system is not 
exactly controllable in H 2 x L 2 (0, 1) with controls in L[ oc ([0, +00)), r > 1. However, it is 
proved in [3] that this system is exactly controllable in H b+ x H 3+ (0, 1) with controls in 
Hq(0,T), at least locally around a stationnary trajectory. Similarly, Turinici's generaliza- 
tion |l33 of Theorem Q] applies to ID Schrodinger equations of the form 



where the state is if>, the control is u and \i € C°°([0, 1]). It proves that this system is not 
exactly controllable in H 2 ((0, 1),C) with controls in L 2 OC ([0, +00)). However it is proved in 
[21 S| that this system, with n(x) = (x — 1/2) is exactly controllable in _ff 7 ((0, 1), C) with 
controls in Hq(0,T), locally around the eigenstates, for T large enough. 

The conclusion of [21 El E] is that, sometimes, the negative result of Theorem Q] is only 
due to a bad choice of functional spaces that do not allow the controllability; but positive 
controllability results may be expected in different functional spaces. Therefore, one may still 
hope to prove the exact controllability of the Bloch equation in some well chosen functional 
spaces. We will see in this article that it is not the case: the Bloch equation is not exactly 
controllable in a much stronger sense than the one of Theorem [U 

Indeed, we will prove that, when (w*, w*) = (—00, +00), the reachable set (in finite time 
and with small controls) from M = is a submanifold of some functional space, that 
does not coincide with one of its tangent spaces. When the domain (w*,u>*) is a bounded 
interval of R, we will see that there exist analytic targets, arbitrarily close to ej, that cannot 
be reached exactly from 63 with bounded controls in L 2 (0,T). Thus, the non controllability 
of © is not related to a regularity problem and this equation corresponds to a very different 
situation from OElll]. 





0, xe (0,1), t £ (0,+oo) 




x e (0,1), te (0, +00), 
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1.3 Outline and open problems 



In section [2 we study the linearized system of ([3]) around the steady-state (M = 63, (u, v) = 
0) with —00 < u>* < ui* < +00. This system is shown to be approximately controllable in 
C°([w*, w*],R 3 ), in any finite time T, with unbounded controls (u,v) € C c °°((0, T), R 2 ). But 
it is not exactly controllable neither in finite time nor in infinite time. Moreover, for any 
reachable target, there exists only one control which steers the control system to the target. 

In section [31 we study the exact controllability of the nonlinear system (J3j) , locally around 
M = €3, in finite time. First, we prove that the simultaneous exact controllability with 
respect to ui in the whole space R (i.e. = —00, ui* — +00) does not hold with bounded 
controls. Indeed, for every time T > 0, the reachable set from M = e 3 with bounded 
controls in L 2 (0, T) is a strict submanifold (of some functional space) that does not coincide 
with one of its tangent space. Then, with an analyticity argument, we deduce that the 
simultaneous exact controllability with respect to u> in a bounded interval (w*,cj*), —00 < 
uj« < ui* < +00, does not hold neither. 

The exact controllability of ((3j being impossible with bounded controls, in sections 0] 
and 03 we investigate the exact controllability of with unbounded controls. 

In section HI completing the arguments of [SUTOl [TT] , we prove the ensemble controllabil- 
ity of (J3j : any measurable initial condition Mq : (oj* ,u>*) — > § 2 can be steered approximately 
in L 2 (lu«, ui*) to 63. This approximate controllability indeed holds for stronger norms, for 
instance and Vs S (0, 1). The controls used to realize this motion are sequences 

of pulses presented in [9] (but one may also use controls in L™ c ([0, +00))) and the proof 
relies on non-commutativity and functional analysis. 

In section [U we propose other explicit unbounded controls realizing the asymptotic local 
(exact) controllability to ea, simultaneously with respect to w in a bounded interval. Here, 
the proof relies on Fourier analysis. 

Finally, in section [6] , we compare the feasibility, the time and the cost of the two 
controllability processes presented in sections H] and on particular motions. 

Let us emphasize that the behavior of the nonlinear system around e% is very different 
from the one of the linearized system around 63. Indeed, 

• first, the linearized system is not asymptotically zero controllable whereas the nonlin- 
ear system is asymptotically locally controllable to e^, 

• then, as seen in section for the linearized system and for any reachable target, only 
a single control works, whereas for the nonlinear system and for any initial condition, 
many controls allow to reach exactly (in infinite time). 

Thus, the nonlinearity allows to recover controllability. Finally, let us mention some open 
problems. 

In section [3j we prove the non exact controllability to e% with bounded controls, in 
finite time, because the reachable set is a submanifold. The equation of this submanifold 
and the validity of the same negative result in infinite time (i.e. the non asymptotic exact 
controllability to with bounded controls) are open problems. 

In section O we prove the exact controllability to e$ with unbounded controls, in infinite 
time. The validity of the same result in finite time is also open. 

Before starting the mathematical study let us introduce some notations that will be used 
in all the paper. We write 




Z(t, uS) := (x + iy)(t, ui), w(t) := (— v + iu)(t). 



(6) 
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Thus, when, for some time T > 0, z(t,uj) > on (0, T) x then the system 

implies 

dZ 

— (t,oj) = iuZ(t,w) - w(t)y/l - |Z(t,w)| 2 , (t,w)e(0,T)x(w*,w*), (7) 

so 

Z(t,w) = (z (u)-J w(T)^/l~\Z(T,cu)\ 2 e-^ T dT^ e™\ (t,w) G (0,T) x 

Unless otherwise specified, the functions considered are complex valued and, for example, 
we write L 2 (M) for L 2 (R 7 C). When the functions considered are real valued we specify it 
and, for example, we write L 2 (R, R). 

2 Linearized system around (M = 63, (u, v) = 0) 

In this section — oo < < u* < +00. We are interested in the linearized system of 
around (M = e$,{u,v) = 0), or, equivalently in the linearized system of Q around 
(Z = 0,w = 0), 

Z(t, lu) = iuZ(t, cj) - w(t), Z(0, u) = Z q {uj), (8) 
whose solution is t 

Z(t,u) = (z a (w) - J w{T)e- iUT dr^e iut . (9) 

We prove its non exact controllability and its approximate controllability with unbounded 
controls. 

2.1 Non exact controllability 

We denote by T the 1-D Fourier transform: 

T{w){uo) = / w{t)e~ lut dt. 
Jr 

When a function is defined on 1 C 1, we extend it by on K \ I. One has the following 
proposition. 

Proposition 1 Let T G (0, +00). The reachable set from Zq — for {3p with controls 
w G L l (Q,T) is 

{Z(T);w&L\0,T)}=^[L\-T,0)]. 

The set of initial conditions Zq that are asymptotically zero controllable with controls 
w G i 1 (0, +00) for {3p is JF[L 1 (0, +00)]. 

For every Zq G T\L (0, +00)], £/ie function w := T~ x \Zq\ is the unique control in 
L 1 (0,+oo) that steers the control system f^) from Zq to 0. 

Proof of Proposition Q} The two first statements are direct consequences of the 
explicit expression J9]). Concerning the third statement, it is sufficient to prove that if 
w G L (0,+oo) and if T[w] = on then w = 0. Let w be such a function and 

consider (p : C + U C_ U (w*, cj*) — > C, defined by 

r fH(z) , if z g c_, 

<p{z) := I T[w](z), if zgC + , 
[ 0, if z G (u>*,cj*), 

where C + := {z G C; 3(z) > 0} and C_ := {z G C; 3(z) < 0}. Then <p is holomorphic on 
C + and on C_ and continuous on C + U C_ U (uj*,uj*), so it is holomorphic on C + U C_ U 
Since ip vanishes on then ip = 0. Thus w = 0. □ 
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2.2 Approximate controllability with unbounded controls 

Proposition 2 LetT > 0, Z f e C°([cj,, u/*]) and rj > 0. There exists w € C c °°((0,T)) suc/i 
i/iai ifte solution of |2p wii/i Zq = satisfies \\Z(T) — ZfW^aa/^^^,-) < r/. 

Proof of Proposition [H Let T > 0, Z f e C°([lj*,uj*]) and 77 > 0. Thanks to the 
Weierstrass theorem, there exists a polynomial P £ C[X] such that 

WZf^e-™ % - P[iu)\\ L «. { „„„.) <|. 

Applying the control u>(t) := —P(dt)6T/2{t) m © with Zo = 0, we get 

Z(T,u) = -^H(o;)e^ T = P{iu)e-^ e lujT = P^wje*"* 

thus, 

\\Z(T) - Zfh^,^ < |. 

Now, let us smooth this control candidate in order to provide a smooth control. Let g € 
<7~((-l,l),R+) such that J R g = 1. For e G (0,T/2), the function 

,(0:=i,(^) 

is supported in (0,T). Applying the control w e (t) := —P(d t )g e {t) in |(8|), we get 

= P(iu>)g e (uj)e iuT . 

Noticing that 

&u>) - = J ^ g(y)[e-^y _ i] d|/ e — S , 

we get 

llPM^M-e-^JII^^^j <c: ||P(*w)|U- ((i; . >llJ . ) ||5e(w)-e- <w T|| ioo(u>i((J . ) _o 
when e — > 0. Thus, for e small enough, 

\\Z(T) - Z/|U~ KjW .) < 

3 Non exact controllability with bounded controls 

In this section, we study the reachable set from M(0) = for (TJ with bounded controls 
£ L 2 ((0,T),R 2 ). Notice that, when M(0) = e 3 and w is small enough in L^O,! 1 ), 
then z(t,u) > for every (t,uj) G (0,T) x R and 

Z(t,u>) = - [ w{r)y/l- \Z(r,uj)\ 2 e- WT dTe luJ \\l(t,uj) e [0, T] x M. (10) 
Jo 

3.1 Case a;* = —00, c<j* = +00 

In this section, we take = —00, uj* — +00. In a first subsection we precise the functional 
framework in which lflO|) is well posed. In a second subsection, we prove that the reachable 
set from zero, with bounded controls (u,v) in L 2 ((0,T),R 2 ) is a strict submanifold of some 
functional space that does not coincide with its tangent space at zero. In particular, ([3]) is 
not locally controllable with bounded controls (u,v) in L 2 ((0, T), R 2 ). 
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3.1.1 Solutions on [0, T] 

Proposition 3 LetT > andR := l/(2\/T). For every w € L 2 (0,T) with |HIl 2 (o,t) < R, 
there exists a unique Z € C°([0, T], L 2 (M)) n C°([0,T] x K) solution of f!0\) and it satisfies 



||^IU°°((0,T)xR) ^ vT||u;||i,2(o,T)5 (11) 



II^IIc°([o,t],l 2 (r)) ^ 2v27r||u>|| i 2( 0j < r ). (12) 
Moreover, for every w\,W2 £ L 2 (0,T) with ||wi||z,2(o,t) < R an d ||^2||z, 2 (o,T) < R> we have 



\\Zi - Z 2 ||z,oo(( ,t)xR) < 2vT||wi - w 2 \\ L ^(a,T), (13) 

\\Z\ - Z 2 \\c°([o,t],l*(r)) < 4V27r||w7i - w 2 \\l2(o,t), (14) 
where, for j <G {1, 2}, Zj denotes the unique solution of (Q7|) for w := wj . 

Proof of Proposition [3l Let T > and c := 1/V3, which is chosen so that 



\f'(x)\ ^c,Vxe [0, 1/2], where /(a;) := yJX-x 1 . (15) 

Let u; G L 2 (0,T) be such that ||«j||i,2(o,t) < R- We apply the Banach fixed point theorem 
to the map 9 defined on 

B := C°([0, T], L 2 (M)) n C°([0, T] x M, Bfc(0, 1/2)) 

by ©(O = Z where 

B c (0,l/2)):=UeC; \z\ < 1/2}, 



Z(t,w) = - / w(t)V1 - |C(r, W )| 2 e-^ T dre Mt ,V(i, W ) G [0,T] x 



Note that B is a nonempty closed subset of the Banach space C°([0, T],L 2 

First step: takes its values in B because R\/T ^ 1/2. 
Let £ € £? and Z := 0(£). The Cauchy-Schwarz inequality leads to 

|Z(t,w)| ^ IMUw) < Vt|MU»(o,T) < VTi? < 1/2. (16) 
Thus Z G C°([0, T] x R, B c (0, 1/2)). Thanks to the decomposition 

Z(i, W ) - -^[r-ttoip,*]]^) - y ^(r)(Vl-|e(T, W )| 2 - l) e-^dre-*, (17) 

where T a ((p)(s) := tp(s-a), thanks to the Plancherel theorem, ((TSj) and the Cauchy-Schwarz 
inequality, we get 



ll^(^)IU 2 (R) ^ V / 27r||w|| £ 2( T) + 

\/27r||w|| L 2( T ) + (j M ||Hlia(0,T) 



tUJT dT 
1/2 



2 \ V2 

dui 



IHU 2 (0,T)(V27T + C\/T||£|| c n([o iT ], L 2( R ))^, 

(18) 

so Z{t) € L 2 (M) for every t e [0,T], In the right-hand side of l[T7|). the first term belongs 
to C°([0, T], L 2 (R)) and the second term also, as one can prove by applying the dominated 
convergence theorem with the following domination, that holds for every (t, u>) € [0, T] x R, 



w(r)^l-\ar,co)\ 2 - l)e— dr| < c|HIl 2 (o,t) / |f(r,w)| 



T \ 1/2 

2 I 
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Second step: is a contraction on B because c^/TR < 1. 
Let 6,6 G B, Zi := 6(6) and Z 2 := 6(6)- We have 

(Z 1 -Z 2 )(t,u) = -J «;(r)(Vl-|Ci(r,w)| 2 - ^1 - ^^^"^re^. 

Using and the Cauchy-Schwarz inequality, we get 

\\Zl - ^2||l~((0 : T)xR) ^ HwlU^O.Tjcll^l — 6IU°°((0,T)xR) ^ cVTR\\^i — 6 II L=° ((0,T) xR) ■ 

Working as in {Ht , we a l so get 

\\(z 1 -z 2 ){t,.)\\ L 2 m <cVri?||6-6llco([o J T],i 2 (R) )l vte [o,r]. 

Third step: Proof of 07p and 02P £/ianfcs fo cv^i? < 1/2. 
Since cVTR ^1/2, the inequalities ifTTj) and (jT2j) are consequences of (fl6|l and lfl8|) with 

.Fouri/i siep: Proof of J23J) and i/ian&s io c^Ti? < 1/2. 
Using the decomposition 

(Zi - Z 2 )(t,w) = - f*( Wl ~ w 2 ){ T )^l-\Z 1 {T ,u J )\ 2 e-^dTe^ t 

-!lMr){^l- \Zx{r,uj)Y- v /l-|Z 2 (r,a;)|2)e-^dre- t , 

the Cauchy-Schwarz inequality and (fT5l) . we get 

|(Zi - Z 2 )(t,w)| < v^||t«i -HU 2 (o,t) +cs/TR\\Z x - Z 2 \\ L ~ mT)xm ,\J(t,uj) G [0,T] x K 
which, since cVTR ^ 1/2, leads to 

11^1 - #2|U°°((0,T)xR) ^ 2\/r||wi - W 2 || L 2 (0,T)- 

Using the decomposition 

(Zi - Z 2 )(t,u) = -T[T- t {wi - w 2 )\m tt ]](cj) 

- fo (wi-w 2 ) (r) ( y/l-lZ^T,*)? - l) e-^ dre iut 

-IoMr)(V l - \Zi(t^)\ 2 - V 1 - |^ 2 (r, W )| 2 )e— dre- 4 , 

and working as in (fl8| . we get 

||(Zi - Z 2 )(t)\\ L 2(^ < \/27r || wi - w 2 ||l2 (0iT ) + ||wi - W2||i2( 0iT )cVT||Zi|j c .o([ 0:T ] :L 2( R)) 

+ \\w 2 \\ L2(0,T)cVT\\Zl - 2 , 2 ||c°([o,t],l 2 (r))- 

Thus, using JUJ) and c\/Ti? < 1/2, we get 

ll^l - Z 2 ||c n ([0,T],L 2 (R)) < 2\/27r^l + 2cVTRj\\Wi - W 2 ||l 2 (0,T) < 4v^7r||wi - W2||i 2 (0,T)- 

This shows the existence the existence part of Proposition [3] and the uniqueness if one 
requires that Z takes its values in Bc(0, 1/2). The uniqueness without assuming this last 
assumption can be easily obtained from the previous study by noticing that this study 
implies that, if two solutions are equal on [0,r] with r G [0, T), then there are equal [0, r'] 
for t' > r. □ 
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3.1.2 Structure of the reachable set from zero in time T 

The goal of this section is the proof of the following result, where 

B R [L 2 (0,T)} := {w £ L 2 (0,T); \\w\\ L 2 {0 ,t) < R}. 
Theorem 2 Let T > and R := l/(4\/3T). The image of the end point map 



F T : B R [L 2 (0,T)} -> L 2 nC»(l) 

it) i ► Z(T, .) where Z solves (f70| 



(19) 

is a sin'ci submanifold of L 2 n C° (R) o/ infinite codimension that does not coincide with its 
tangent space at zero. 

The proof of Theorem [2] relies on the following results (see [UJ Theorem 73. E and 
Corollary 73.45, Chapter 73]). 

Theorem 3 Let M and N be two C k -Banach manifolds with chart space over R and k £ 
{oo} U N \ {0}. Let F : M ^ N be a map of class C k . If F is a C k embedding, then 
S := F(M) is a submanifold of N and in particular a C k -Banach manifold. 

Theorem 4 Under the same assumptions as in Theorem^ if F is an injective C k immer- 
sion and if F is closed, then F is a C k emb 



Proof of Theorem [2} We take M := B R [L 2 (0,T)} and N := L 2 n C£(R). They are 
both C°°-Banach manifolds as open subsets of Banach spaces. The continuity of Ft is a 
consequence of (fl3| and lfl4|) . With similar manipulations as in the proof of (fl3| and l|14p , 
one can prove that Ft is C 1 and dFT{w).W = £(T, .) where £ is defined by 



f (t, w) = - /„ W(t) y /1 - \Z( r,uj)\ 2 e- luJT dTe iuJt 

+ Jo* «>(t) * [ f (T ' " )€(r ' " )] e-*W* V(t, a,) e [Q, T] x R. (20) 
^l-|Z(r,w)| 2 

We use the same notation c as in the previous proof (see (flSj) ), 



Krsi siep: Ft «s injective on B R [L 2 (0, T)] because Qc^/TFL < 1. 
Let w 1 ,w 2 £ B R [L 2 (0,T)} be such that F T (wi) = F T {w 2 ). From 



T 

w l {t)y/l^\Z^t^e- iut dt= I ^(tjVl-l^^wJPe-*-**, 



'0 

we deduce 



T[ Wi -w 2 ](lu) = J^(w 2 - wx){t)(^/l -\Z 2 (t,Lu)\i- ije-^dt 
+ S*wi(t){ } /l-\Z 2 {t,u)\*- jl-\Z^t,u)\< 

Considering the L 2 (R)-norm of both sides, using Plancherel equality and working as in (jT8 
we get 



27r||wi - W 2 \\l*(0,T) < \\w 2 - Wi\\ L 2( 0> t)CVT\\Z 2 \\c°([0,T],L2(K) 
+ ll w l|U 2 (0,T)C\/T||2'l - Z 2 \\ C 0([ Qt T],L 2 ( 

Using |fl2|) and (fl4"|). we deduce 

\/27r||wi - w 2 ||l 2 (Ut) < 6cV27r\/Ti?||w;i - U^IU^o.t)) 
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which gives the conclusion, because Qc\/TR < 1. 

Second step: Ft is an immersion because 6cVTR < 1. 
Let w G B R [L 2 (0,T)] and W G L 2 (0,T) be such that dF r (io).W = 0. Thanks to JUJ), we 
have 

F[W]{")= ~ So W(t) ( yi - \Z (r, - l) e-^dr 

Jo ^l-\Z(r,u;)\ 2 
Considering the L 2 (R)-norm of both sides and working as in lfT8|) we get 

V2tt\\W\\ L 2^ 0T ) ^ ||W|| i 2 (0 ^ T) cVr||^||c ([0,T],L 2 (R)) + IHU 2 (o,t) c ^II£IIc°([o,t],l 2 (r)) 
Admitting the following inequality, 

I|£||c°([o,t],l 2 (r)) < 4\/27r||W r || i 2(o i T), (21) 

and using lfT2|) . we get 

L 2 (0,T), 



which gives the conclusion because 6c\/TR < 1. Now, let us prove l|21j) . Using the decom- 
position 

£(t,w) = -J-[r_ t W]( W ) 

-So W{s)Lfl - \Z{s,uj)\ 2 - lje-^dse^ 



Jo y/l-\Z( S ,u;)\* 

and working as in ((18)1 . we get 

||£(£)IU 2 (r) ^ V27t||W / ||l 2 (o : t) + ||W||l 2 (o : t)cVT'||^||c ([o,t],l 2 (r)) 

+ IMU 2 (0,T)<VT||£|| c O([ 0: t],L 2 (R))- 

Using mi) and cVTR < 1/2, we get {H}. 

Third siep: i*r is a closed map because 6cVTR ^1/2. 
Let A be a closed subset of -B^[L 2 (0, T)]. Let (Z n (T,.) = i 7 r(u' n )) ne N be a sequence of 
F T (A) that converges to ^oo(-) in L 2 n C°(R). In order to prove that € -Ft (A), we 
prove that (w n ) n ^ is a Cauchy sequence in L 2 (0,T). For every n £ N, we have 

Z n (T,u>) = -T[T_ T w n ]{u) ~ J w n {t)(^fl-\Z n {t,u)\ 2 - iy-^dte^ T , 

so, for n,p£ N, we have 

^[r-T^n - w p )](uj) = (Z p - Z n )(T,u) 

- / T K - Wp ) (t) ( v/T= |Z„(^)| 2 - l) e-^dte^ 

- / T (yi-|Z„(t, w )| 2 - Vl-li^w)! 2 ) e-^dte*- 2 ". 

Considering the L 2 (R)-norm of each side, using the Plancherel equality and working as in 
(HI), we get 



2i:\\w n - W p \\ L i(p,T) < IK^n ~ ^)( T )IIl 2 (R) 

+ ||w„ - w p ||l 2 (o,t)cVT||Z„|| c o([ 0! t] : l 2 (r)) 

+ ll w plli 2 (0,T)CV / T'||2'ri - ^p||c°([0 : T],L 2 (R))- 
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Using (TU and ([Tijl. we get 

V27r||w„ - w p || L 2 (0:T ) ^ \\(Z„ - Z p )(T)|| L 2 (R) + 6c\/27r\/'Z ; R||w, i - w p || i 2 (0:T) . 

which gives the conclusion because Qc^/TR =1/2. 

Fourth step: The manifold S :— FtBr[L 2 (0,T)] does not coincide with its tangent space 
at zero. 
We have 

dF T (0).W = -T[t_ t W\ 

thus 

T S = T[L 2 ({-T,0))]. 
Let us compute the third order development of Ft around 0, 

w (t) = eWi(t) + e 2 W 2 (t) + e 3 W 3 (t) + ... 

Z(t,u) = eZi(t,u) + e 2 Z 2 {t,u;)+e 3 Z 3 (t,uj) + ... 

Since, as e — > 0, 



^l-\eZ 1 + e 2 Z 2 + e^Z 3 \ 2 = ^1 - e 2 \Z x \ 2 + o(e 2 ) = 1 - -e 2 \Z,\ 2 + o(e 2 ), 

we have 

Z 1 (t,w) = -^[r_t(Wi)| [0 , t] ](w), 
Z 2 (t,uj) =-JF[ T _ t (W 2 )| [ o > t]]M, 



Z 3 (t,o;) = -^r-t(Wa)| [0l t]]H - 5 / ^(r^^a^e"^^^ 4 



We want to prove the existence of Wi S L 2 (0,T) such that the map 

i-T 

2 ■— iuir , 



W 1 (T)|Zi(r,«)|^e- ,WT dT 

does not belong to .F[L 2 (0,T)]. Using the explicit expression of Z\, the change of variable 
c 2 — > x = t -\- o~x — o~ 2 and Fubini's theorem, we get 

= So W i( T )( So ^(^^^^(lo W 1 (a 2 )e i ^da 2 y- iUT dT 
= £=0 SS 1= o XJ 2 =o ^ 1 (r)lU 1 (a 1 )W 1 Me--(^-^)da 2 da 1 dr 

r T f T rcri+T , 



where 



/r=o Ci T 4o,x-. } WtfWfa) W(r + oi - x)d CTl dr if * e (0, 2T) 
if x g (0,2T). 



Computing the map $1 associated to W = 1[o,t], we get $i(3T/2) = T 2 /16 thus, for e small 
enough, F T (elp )T ]) ^ T S. □ 
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3.2 Case — oo < lu* < u* < +00 

The goal of this section is the proof of the following result. 

Theorem 5 (i) Let T > 0, u,v S L 2 (0,T) and M be the solution of 

^f(t,uj) = [ujQ z + u(t)n x + v(t)n y }M(t,ui), (t,w)e(0,T)xC, (22) 



M(0,cj)=e 3 . 

Then wgCh> Z(T,w) is holomorphic. 

(ii) Let T > and i? := l/(4v3T), There exists Zf : (a;*,w*) — + C analytic such that, 
for every e* > 0, i/iere era'sis e € (0,e*) such that, for every w € -Br[L 2 (0, T)] ; ifte solution 
of (pj) satisfies Z(T) ^ eZ/. 

As a consequence, there are arbitrarily small analytic targets on (a;*, uo*) that cannot be 
reached exactly in finite time, with controls having a prescribed L 2 -bound. 

Proof of Theorem [H (i). Let T > 0, u,v e L 2 (0,T) and M be the solution of <j22j) . 
We introduce the functions M U M 2 :[0,r]xlxl^l 3 defined by 

Mi^wx.wa) :=$R[M(i,wi+iw 2 )], M 2 (t,LJ 1 ,cj 2 ) := R[M(t, wi + iu*)]. 
The function M(t, u>i, uj 2 ) ■= (Mi(t, wi, u> 2 ), M 2 (t, uii, cj 2 ))* solves an equation of the form 

^ = /(t,M ) « 1 ,w 2 ). (23) 

The map / is measurable, of class C 1 with respect to (M, wi, w 2 ) £ M 6 x R x R and satisfies 

|/(t,M,wi,W2)| < (|wi| + |u a | + |u| + 
1/^(4, M,wi,w 2 )| < |wi| + M + |u| + |v| 
\f Ul (t,M, ^,^2)] < |M|, !/„,(*, M,wi,W2)| < |M|. 

Thus M has partial derivatives with respect to oj\ and uo 2 . (To check that, one can, for 
example, adapt the proof of [8j Theorem 4.1, chap. 4, page 100].) Let us prove that they 
satisfy the Cauchy-Riemann relations, in order to get the holomorphy of u> € C 1— * M(T, u>). 
We introduce the notation 

r fc ,i (*, wi , wa) := ^ (t , wi , ua) , for k,l £ {1, 2}. 

Differentiating the system f23|) with respect to oji and w 2 , we get 

d ^~ Y ^ = A(t,ui){Y u - y 22 ) - B(u, 2 )(Y 12 + Y 21 ), 
0(y 12 +y 21 ) =A(t,Lj 1 )(Y 12 +Y 21 ) + B(u 2 )(Y 11 -Y 22 ), 
{Y 11 -Y 22 ){0,lu 1 ,lu 2 ) = 0, 
(*ia+l r 2i)(0,wi,w 2 ) = 0, 

where 

-wi v(t) \ / -w 2 

A(t,u;i):=| wi -u(t) and B(w 2 ) := lu 2 

-v(t) u{t) / \ 

The uniqueness of the solution of the Cauchy problem ensures that In = F 22 and 
Y\ 2 = —Y 2 \. 
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(ii) Let Zj : R — > C be an analytic function that does not belong to the tangent space 
of the image of Ft at zero (i.e. which is not the Fourier transform of an a function in 
L 2 {{— T, 0))). Then, for every e* > 0, there exists e G (0, e*) such that eZf does not belong 
to the image of Ft- Thanks to (i), reaching eZf on (a;*, to*) in time T with controls u and 
v in L 2 ((0,T),R) is equivalent to reaching it on R. But eZf does not belong to the image 
of Ft- Therefore eZf cannot be reached on (a;*, a;*), in time T, with controls u and v in 
B fl [L 2 ((0,T),R)]. □ 

4 Ensemble controllability with unbounded controls 

In this section, we take — oo < uj* < u* < +oo. The goal of this section is to complete 
the very interesting arguments of [9l [TUl [H] with functional analysis ideas, to prove the 
ensemble controllability of ((3]), (i.e. the approximate controllability of |(3]) in L 2 (uj*,uj*)) 
with unbounded controls. Actually, we prove a stronger result. 

First, let us introduce the definition of solutions of |(3|) with Dirac controls. 

Definition 1 Let b G [0, +oo), /3, 7 G R and M : (ui*,lu*) — > § 2 . The solution of (TJj with 
M(Q) = M , u(t) = f3S b (t), v(t) = j8 b (t) is 



M(t,Lu) = 



exp(ujn z t)M (uj) for t G [0,6), 

exp(u)Q z (t — b)) exp(f3Q x + jQy) exp(u)fl z b)Mo(ui) for t 6 (6, +00), 



M(b + ,uj) = exp(pn x +-yn v )M(b-,u). 

Let 

H 1 ^,,^*)^ 2 ) := {M G w*),R 3 );M(u) G § 2 ,Vw G (w*,w*)}. 

Let also C/[i; w, u, M ] denote the value at time t, of the solution of ([3]), with initial condition 
Mq at time 0. Thus, U[t; u, v, M$\ is a function of u> G (w*,o;*). Definition [T] is motivated 
by the following result, which is a consequence of explicit expressions and the boundedness 
of (cj*,u;*). 



Proposition 4 For every (3, 7 G K, we /lave 



lim 



£7 



/3 7 

^ + e ! — l[b,6+e]i ~l[b,6+e 



U[b+-I38 bll 8 b ,. 



= 0. 

£(JT 1 ((a)»,a)*),R 3 ),JT 1 ((w,,a)»),R 8 )) 



Let us introduce the set D of finite sums of Dirac masses on [0, +00). The goal of this 
section is to prove the following result. 

Theorem 6 Lei M G H 1 ^*, w*), § 2 ). There exist (t n ) neN G [0, +oo) N , and 
(u n )neN, (vn)neN G D N such that 

U[t^', u n ,v n , Mq] — > e3 weakly in H 1 ((w*, u;*), R 3 ). 

Thanks to Proposition [4j one easily get, from Theorem [6] the following corollary. 

Corollary 1 Let M G if 1 ^*, w*),S 2 ). T/iere ea;isi (* n )neN e [0, +oo) N , and 
K)n£N, (fn)neN G ( [0, +00) , R) N smc/i £/ia£ 

U[t M ] — > e 3 weakly in H ((w*, o/),R 3 ). 
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Thanks to the compactness of the injection H 1 ^*, to*) — ► L 2 (lu*,uj*), and the time 
reversibility of ([3]), Theorem [6] and Corollary Q] give the ensemble controllability of 
(i.e. its approximate controllability, for the L 2 ((w* , to*), M 3 )-norm, in finite time). These 
statements also give the approximate controllability of J3]), for the L°°((o;*, w*),R 3 )-norm, 
in finite time. The proof of Theorem [6] relies on the following Lemma, that will be proved 
later on. 



Lemma 1(1) Let M £ H 1 ^, to*), S 2 ) be such that M' ^ 0. There exist T > 0, u, v £ D 

• one /ias 

||l/[T + ;«,t; > Af] / || i;a < \\M'\\ L 2, 

• for every sequence (M n )„ e N £ H 1 ((ui* , to*) ,§ 2 ) N satisfying 

\\M^\\ L2 ^\\M'\\ L2 ,Vn£n (24) 

and 

M n -> M wea% in ff^w,, w*),R 3 ) (25) 
i/iere erriste an extraction ip such that 

\\U[T+-u,v,M v(n) r\\ L 2 < ||M; (n) || i2 ,VneN. 

fi-Jj Lei M S § 2 6e shc/i i/iai M ^ 63. There exists £ [0,27r) such that for some 
(u,v) £ {(ttJi + (7T + 9)6 2 ,Q), (0,7ro\ + (tt + 6»)<5 2 )}, £7[2+; «, v, M] is constant over (to*, to*) 
and 

\U[2 + ;u,v,M]-e 3 \ < \M - e 3 \. 

In section |4~T| we prove Theorem [6] thanks to functional analysis and Lemma [IJ which is 
proved in section l4~2l and [4751 In section l4~2"t we recall a preliminary result, which is already 
presented in [9j [TT]. In section [473} we deduce the proof of Lemma [TJ 



4.1 Proof of Theorem [6] thanks to Lemma Q] 

In this section, we deduce Theorem [6] from Lemma [IJ using similar arguments as in [13] . 

Proof of Theorem \6\ thanks to Lemma[TJ Let Mo £ i? 1 ((w >t ,[jj*),S 2 ) be such that 
Mo ^ e 3 (otherwise t n = gives the conclusion). We introduce the set 

K:= {M£H 1 ((to*,to*),S 2 ); 3(i„)„ eN £ [0, oof, 3(u n ) neN , (v n ) neN £ D N 
such that \\U[t+;u n ,v n ,Mo}'\\ L 2 < \\M' Q \\ L 2,Vn £ N 
and U[t+;u n ,v„,M ] -> M weakly in R 3 )} 

and the quantity 

m := inf{|| Af'|| L 2; M G if}. 
Notice that if is not empty because it contains M$ (take t„ = 0). 

First step: Let us prove the existence of e £ K such that ||e'||x,2 = to. 

Let (M n ) n6 N* £ K n be such that IjM/Jj^ — > to when n — > +00. Then (M„) n£ N* is a 
bounded sequence in i? 1 ((w*,u;*),K 3 ), thus, there exists e € if 1 ((a;*, a;*), § 2 ) such that (up 
to an extraction) 

M n — > e weakly in ii 1 and strongly in L 2 . (26) 

Then, 

||e'|| L a <Uminf||M;|| i2 =m. 
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Let us prove that e belongs to K, which gives the conclusion of the first step. For every 
neW, M n e K so there exist (^) pG n € [0, +oo) N , (uP) peN , {v^) p&i € D N such that 

\\U[tP n + ;uP,vP,M ]'\\ L 2 < HM^.VpGN.VnGN* (27) 

and 

u^, Mo] — > M„ weakly in if 1 and strongly in L 2 , when p — > +00, Vn e N*. 
For every n e M*, we choose p = € N such that 

||^W+ ;u P(n) ^(n) M()] _ Mn || ^ I. (28) 

n 

The sequence (y„ := J7[t^ j M)])neN* is bounded in i? 1 ((ti'*,CL'*),M 3 ) because 
of lf27j) . Thus, there exists eg S if ((a;*,a;*),§ ) such that (up to an extraction) 

Y n — ► ej weakly in ii 1 and strongly in L 2 . 

The definition of K ensures that ej belongs to K. Moreover, because of (f26| and l|28jl . 
Y n — » e strongly in L 2 ((w*, K 3 ), thus (uniqueness of the strong L 2 limit), e = ej and 
e S K . 

Second step: Let us prove that m = 0. 

Working by contradiction, we assume that m > 0. Then e' ^ thus, we can apply 
Lemma [T](l). There exist T > 0, u,v £ D such that 

\\U[T+;u,v,e}'\\ L 2 < \\e'\\ L 2 = m, 

and an extraction tp such that, with the notations of the first step, 

\\U[T+;u,v,Y v(n) ]'\\ L2 ^ ||r; (n) || £a ,Vn G N. (29) 

Let us prove that U[T + ;u,v,e) belongs to K, which gives the contradiction. 
Using and ([2"9]t. we have 

||f/[T+; W , U) F v , ( „ ) ] , || i2 < ||M || L2 ,VneN. 

Thus, there exists 6 iJ 1 ((w >t ,cj*),S 2 ) such that (up to an extraction) 

U[T + ; u, v, 1^(„)] — > e* weakly in ii 1 and strongly in L 2 . 

Then, e* G if, by definition of if. But we have 

\\U[T + ;u,v,Y v ( n )} - U[T + ;u,v,e]\\ L 2 = \\Y^ n ) - e|| L 2 — > when n -> +00, 

thus £/[T + , it, v, e] = e* (uniqueness of the strong L 2 limit), and U[T + , u, v, e] € if. This 
ends the proof of the second step. 

Third step. With a slight abuse of notations, let us still denotes by § 2 the set of constant 
functions from (— lu*,lu*) with values into § 2 . Thanks to the first and second steps, the set 
K n § 2 is not empty, so we can consider 

m := inf{|M- e 3 |;M e ifnS 2 }. 

Working exactly as in the first step, one can prove that K nS 2 is a closed subset of S 2 , thus 
K n § 2 is compact and there exists eeifflS 2 such that |e — e3 1 = fh. 
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Fourth step: Let us prove that m = 0, which gives the conclusion. 

Working by contradiction, we assume in > 0. Then e ^ e% and we can apply Lemma Q] 
(2). There exists 9 G [0, 2vr) such that, for some (u,v) G {(ttSi + (tt + 6)5 2 , 0), (0,tt£i + (tt + 
8)62)}, U[2 + ;u,v,e] is constant over (u>*,ui*) and 

|f7[2+;u,w,e]-e 3 | < |e-e 3 |. 

Let us prove that U[2 + ;u,v, e] belongs to K, which gives the contradiction. 
First, we emphasize that explicit computations show that, 

exp(7rl7^) exp(wf2 z ) exp(7ril^) = exp(— u>Q, z ), V£ G {x,y}. (30) 

Thus, for some £ G {a;,y}, we have 

U[2 + ; w, u, .] = exp((# + 7r)i\) exp(uiil z ) exp(7rf2^) exp(ujil z ) 

= cxp(9fl^) exp(nQ^) exp(wf2 2 ) exp(nQ^) exp(cjf2 z ) 
= exp(8il^) exp(— luQ z ) exp(ujtt z ) 
— exp(6T^). 

Since e £ K, there exist (s n )„ G N € [0, +oo) N , (/in)neN, (^)nes G -D N such that 

||C/[ S +;/i n ,^,Mo]'|| L2 < ||A^|| za ,VneN, 

J7[s+; jti n , f n , M ] — > e weakly in H 1 . 
Let Z„ := J7[s+;/i„,^ n ,Mo]. We have 

U[2 + ;u,v,Z n ]=exp(9Slt)Z n . 

Thus 

||C/[2+; U ,«,Z„]'|| i2 = \\Z'J L 2 < ||A/£|| L 2,VneN 
and U[2 + ;u,v,Z n ] -> exp(0fi € )e = £/[2+; u, v, e] weakly in H l . Thus E/[2+;u,u,e] G 

4.2 The argument of P flT] 

The goal of this section is to recall the proof of the following result, which is already presented 
in [21 [□]. 

Proposition 5 Let P,Q& M[X]. The flow of (E) can generate 

I + r{P(uj)tt x + Q(tj)Qy] + o(r) when r -> 0, 

controls that are finite sums of Dirac masses. More precisely, for every e > 0, i/iere 

exists t* = T*(P,Q,e) > such that, for every r G [0,r*], there exist T > and u,v € D, 
such that 

\U[T+; u, v, .]-(/ + t[P(w)Q* + Q(w)n y ]) 

Remark 1 Let us explain why Proposition [H may not be sufficient to prove the ensemble 
controllability (i.e. the approximate L 2 (ui^,bj*)- controllability) of with Dirac controls. 

First, let us remark that, for every point M — (x^\ x^ 2 \ x^) G § 2 , such that =/= 0, 
then (fl x M,rt y M) is a basis ofT§?M (the tangent space o/S at M). 

Let e > and Mq = (xo, yoi z o) G L ((a?*, u>*), § 2 ) be such that, z(u>) 7^ for almost 
every w G (u>*,ui*). Following a classical strategy, we consider an homotopy 

H : [0,1] x -> § 2 

(s , w) 1— * H(s,u>) 



«C er. 
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such that H G C^QO, 1], L 2 ((w*, uj*), S 2 )), H{0,u) = M (lu), H(1,w) = e 3 , and we try to 
reach e% from Mq by following the path given by H. Since z ^ a.e. on (u)*,ui*), there exist 
f,g€ L 2 ((lu*,uj*),R) such that 

OH 

— (0,w) = f(u)Q x M (u) + g(u;)QyMo(u). 
os 

Thanks to the Weierstrass theorem, there exist P,Q £ WL[X] such that \\f — P\\l 2 < £ and 
\\d — Q\\l 2 < e - Applying Proposition^ one may follow (approximately) the direction given 
by Q^-(0,uj), with a small amplitude r* , that depends on this direction. 

If one wants to be sure to reach e% in finite time, by iteration of this process, one would 
need, at least, the independence of the amplitude r* with respect to the direction ( otherwise, 
one may stop in the middle of the path). However the maximum amplitude t* given by 
Proposition^ depends on the direction, through the choice of the polynomials P, Q. 

Proof of Proposition Q>} In this proof r is a positive real number. By (|30|). 

exp(ir£l x ) exp(wf2 z r) exp(— 7rO x ) = exp(— w$7 z r), (31) 

and this evolution is generated, in time r, by the controls u{t) — —■k6q{€) + irS T (t), v = 0. 
The controls 

u(t) = V^oW - (it + -fi) + ^a^f (*) 
(resp. u(t) = -Tr5 {t) + (tt - y/^S^t) + y/rS^t)), 
v = 0, generate in time 2-y/r the evolution 

Ui(t) := exp(irQ x ) exp(u)Q, z y/r) exp(— (tt + \/t)Q x ) exp(u)Cl z ^/r) exp(fl x y/T) 
= exp(-ujfl zy /r) exp(-fl xy /r) exp(wfi 2 i/F) exp(f2 xV / r) 
= I + tuj\Q, z ,Q. x ]+o(t) 

= I + TUjfly + o(t) 

(resp. 

Ux(-t) := exp(n xy /r) exp(ujtt z ^/r) exp((7r - ^/t)Q x ) exp(uj£l z ^) exp(-7rfi x ) 
= exp(Q x y/r) exp(wf2 2 ^r) exp(-tt x ^/r) exp(-ujQ zX /r) 

= I - TLLt[Q z ,tt x ] + o(t) 

= I — TLofly + o(r)) 

where we have used ((3Tj) to pass from the first to the second line. Here and in the following, 
o(t) denote quantities which tend to in the C(H l ((u> x , uj*), R 3 ), i/ 1 ((w*, uj*), R 3 ))-norm as 
t — > + . In the same way, there exist controls, that are sums of Dirac masses, that generate 
in time the evolutions 

U 2 (t) := exp(-wf2 z ^/r)J7i(-r) exp(wrj 2% /T)C/i(r) 
= I + T 3 / 2 uj 2 [n y ,n z }+o(T^ 2 ) 

= I -T 3 / 2 LU 2 n x +0{T 3 / 2 ), 

U 2 (-t) := U\{t) exp(wfl 2 r )[/i(~r) exp{—u)Q. z y/r) 

= i-T 3 / 2 uj 2 [n y ,n z } + o(T^ 2 ) 

= I + T S / 2 LJ 2 n x +0(T 3 / 2 ), 

and in time 14-y/r the evolutions 

f/ 3 (r) := exp(-wil 2 ^/r)J72(-T)exp(w^ 2 \/T)C/2(T) 
= /-r 2 c; 3 [f7 2 ,r2 2; ]+o(r 2 ) 
= /-r 2 c; 3 ^ + o(T 2 ), 
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U 3 (-t) := U 2 {t) exp(wO z r )[/ 2 (-r) exp(-w^ 2 ^/r) 
= I + T 2 LU 3 [fl z ,n x ] + o(t 2 ) 
= / + r 2 o. 3 r! y + o(r 2 ), 

Thus, one can generate I ± tuj 2 £I x + o(r) in time 6r 1//3 , and / ± rw 3 i7 a + o(r) in time 
14-r 1 / 4 . Iterating this process, for every n£N, one can generate /iro; 2 "^ +o(r) and I± 
rcj 2ll+1 r2j / + o(r) in a time T„ that behaves like 4"t^\ The same argument with Cl x replaced 
by fly in Uj(r), j ^ 1, shows that for every n € N, one can generate / ± TLu 2n+1 fl x + o(r) 
and / ± TL0 2n u! y + o(t) in a time T n that behaves like 4"tA: . Thus, for every P, Q € 
one can generate 

/ + t[P(u)Q x + Q(Lu)Q y ] + o(r) 
in finite time, by composing the previous evolutions. □ 

4.3 Proof of Lemma Q] 

The goal of this subsection is the proof of Lemma [U thanks to the previous subsection. 
Proof of Lemma fit 

Proof of (1) of LemmaU It is sufficient to prove this statement under the additional 
assumption 

z^O. (32) 

Indeed, let us assume that it is proved when 11321) holds. Let M — (x, y, z) € iP((a;*, uj*), § 2 ) 
be such that M' ^ and z = 0. Then x 2 + y 2 = 1 thus x ^ or y ^ 0. Let us assume, for 
example, that y ^ 0. Thanks to l[30|). we have 



U[2; ^-5 + ndi, 0, .] = exp(w£7 z ) exp^f^) exp(uSl z ) exp^O^) 

= cxp(w£! z ) exp(irQ x ) exp(cuQ z ) exp(irQ x ) exp(-|J7 :z 
= cxp(cjfi 2 ) exp(— £jf2 z ) exp^f^) 
= expd^). 



Thus the function 
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U[2; — Sq+ttS^M] 




has a non vanishing third component and the L 2 norm of its derivative is the same one as 
M . Applying LemmaU] (1) to U[2; ^-So + irSi , 0, M] , we get the conclusion of LemmaU] (1) 
for M. 

Let M = (x,y,z) S H 1 ^, w*),§ 2 ) be such that M' ^ and z ^ 0. 

First step: Let us prove the existence of P,Q G a > and Tq > shc/i i/iai, /or 

euen/ r <E (0,t *), 



• one /ias 



d r 



/ + r[P(w)fi a: + QHfiy] ) M 



s$ \\M' 



rat, 



(33) 



• for every sequence (M„) ii£ n G if 1 ((a;*, cj*), § 2 ) n satisfying (2~4\ ) and $25)) . there exists 
an extraction tp such that 



d 



duj 



I + T[P(u))fl x + Q(ui)Cl y ] ) M ( 



L 2 



^ ||Af' (n) ||| 2 -™,VneN. (34) 
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Developing the square, we get, for r > and P,Q £ 



du> 



I + r[P(tu)n x + Q(u)Cl y ] ) M 



dui 



I + T[P{u)n x + Q(w)Sly] )M, 



L 2 



L 2 



M'\\ 2 L2 + 2rA(P, Q) + t 2 B(P, Q), 



\M^\\ 2 L2 +2TA n (P,Q)+T 2 B n (P,Q), 



where A(P,Q), A n (P,Q), B(P,Q), B n (P,Q) are real constants. Straightforward computa- 
tions give 

A(P, Q) = £l (£ [[P(u)n x + Q(w)Qy]M{u)] , M'(w))du 
= /£" (p'[-zy' + yz 1 } + Q'[zx' - xztydu. 

We look for P, Q e R[X] such that A(P, Q) < 0. Since A is a linear form in (P, Q) it is 
sufficient to prove that A^ 0. Working by contradiction, we assume A — 0. Thanks to the 
density of polynomials in I/ 2 (w*, u>*), we have 



zy — yz' — 0, zx — xz = 0. 



(35) 



Let / be a nonempty connected component of {u> S (a?*, a;*); z(w) ^ 0}. Since z / 0, such 
a J exists. By ((35]) , there exist a,6eR such that = az(ui) and y(w) = &z(w), Voj S P 
Since M takes values in § 2 , we have 

1 = x(oj) 2 + y(uj) 2 + z(ojf = (a 2 + b 2 + 1)z{uj) 2 . 

This shows that I = (oj*,uj*) and that M is constant over (w*,oj*), which is in contradiction 
with the assumption M' ^ 0. Therefore, there exist P, Q S R[X] such that A(P, Q) < 0. 

For every n G N, we have 

A„(P, Q) = P'i-ZnVn + VnZ'n) + Q 'W ~ ^n^n)- 

Thanks to ((25]), there exists an extraction y> such that 

M v ( n \ — > M weakly in P 1 and strongly in L 2 . 
Then ^4 ¥ ,( n )(P, Q) — > A(P, Q) when n — > +00. Thus, we can assume that 

^(»)( P ' ( 3) < \A{P,Q)yn e N (36) 
(otherwise take another extraction). We have 



y/B n (P,Q) 



£[p(tu)n x + Q(iu)n y ]M 



L- 



P'Wl 2 + UQ'IU' + [||P||l- + ||Q|| L »]||M'|| L 2, 



£[p(cu)n x + Q(tu)n y }M ri 



L 2 



p'\\l> + \\Q'\\l* + [\\p\\l~ + \\Q\\l~]\\mu\ 



L 2 

\p'\\l 2 + WQ'h 2 + [\\p\\l- + IIQIMHM'IU*. 

Let r * = r *(M) > be such that 

to (\\p'\W + WQ'h 2 + [\\p\\l°° + \\q\\loo}\\m'\\ l ,) 2 < |A( ^' Q)I . 



Then, for every r e [0, r *], we have (01 and {34]) with a := -A(P, Q). 
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Second step: Conclusion. 

Let P, Q be as in the first step. Let ei > be such that 



ei\\M\\ H i < 



a 



2\\M'\\ L 2 



(37) 



Let t* — T*(P,Q,ei) be as in Proposition [5] and r* := minjr*,^}. Thanks to Proposition 
[3 there exist T > 0, u, v € L~ c ([0, +oo), R) such that 



[7[T+; u, v, .]-(/ + r* [Pfw)^ + QMfi„]) 
Then, using (j3Sj) , 02J and J37J), we get 



£(ff 1 ,ff 1 ) 



eir x . 



(38) 



||tf[T+;u,u,Af] / || £J i £[^[T;«,w ) Af]- (z + TtPHfla, +Q(w)n s ])M 



do; 



/ + r[P(w)0 K + Q(o;)O a ] M 



< e 1 r 1 *||M|| Jf i + (\\M 



L 2 - "'1 



1/2 



L 2 



L 2 



< \\M'\\&. 



Similarly, we have 



'p(-n) Hi 2 



\\U[T+;u,v,M v(n) ]'\\ L * < ei7?||M v(n) || H1 + ||A/; (n) || L 2 - 5^ 

< £l7f ||M|| H x + ||M; (n) || La - 51I ^' 

<n^; ( n)iu- 

This ends the proof of the first statement of Lemma [TJ 

Proof of (2) of Lemma [H Let M = (x, y, z) £ § 2 be such that M ^ e 3 . Then x ^ 
or y 7^ 0. We assume, for example, that y ^ 0. Thanks to f30|) . we have 



U[2 + ; ttSi + (tt + l9)J 2 , 0, .] = exp((-7r + 9)Q X ) exp(ujil z ) exp(irQ x ) exp(uf2 z ) 

= exp(9£l x ) exp(irfl x ) exp(a>f2 z ) exp^f^) exp(u>il z ) 
— exp(0£l x ) exp(— u>Q, z ) exp(wf2 z ) 
= exp(0fi x ). 



Thus, 



1 

U[2 + ; ndi + (tt + 0)5 2 ,0, M] = I cos(<9) -sin((9) I M. 



sin(0) cos(0) 



We get the conclusion with 8 £ [0, 2n) such that 

C/[2+;7r<5i + (tt + 6>)c5 2 , 0, M] = 



x 




.□ 



5 Explicit controls for the asymptotic exact controllabil- 
ity to e 3 

In this section, uj* = 0, uj* = tt. We propose explicit controls realizing the asymptotic exact 
controllability to — es, locally around —e^. 
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First, let us introduce some notations. For a function / : (— 7r,7r) — > C, we denote by 
c n (/) its Fourier coefficients and by N(f) their ^-norm: 



H^e-^dw, N{f) := £ Mfll- 



For a function / : (0, n) — > C, we define 

c n (/) := c„(/),Vn e Z, iV(/) := N(f), 

where / : (— n, n) — ► C, f(u>) ■= f(\u>\). For a vector valued map M = (x, y, z) : (0, tt) — > R 3 , 
we define N(M) := N(x) + N(y) + N(z). Then, we have the following results. 

Lemma 2 For every /, <? : [0, 7r] — > C smc/i i/iai N(f), N(g) < oo, we ftave 

^(/.9) < N(f)N(g). (39) 
For every (x, y) € ^ 1 ((0, n), R 2 ), we /lave, with Z := x + iy, 

i (iV(x) + 7V(y)) < 7V(Z) < iV(x) + N(y). (40) 

For every M : [0,7r] — > S 2 smc/i i/iai N(M) < +oo and z(w) > 0,Vw S [— 7T, tt], we /lave 

iV(z - 1) < 2N(Z) 2 . (41) 

//, moreover, N(Z) ^1/4, i/ien 

\n(Z) < JV(M - e 3 ) < 3iV(Z). (42) 

As a consequence, for a map M : [0,7r] — > § 2 such that N(Z) $C 1/4 and z > 0, the 
quantity N(Z) measures the TV-distance from M to e3. 

Proof of Lemma [2} We have 

N (fa) = E i E c -p(/) c p(ff)i < E M^)i E i^-ptf)! = ^(/Ms)- 

The inequality (f40| is a consequence of the triangular inequality because N(Z) — N(x + iy) 
and iV(x) + N{y) = N((Z + Z)/2) + N((Z - Z)/2i), 

Let M : [0, tt] -> § 2 be such that N(M) < +oo and z > 0. We have 

j\T(z - 1) = -L 1 - ^1-|^MI 2 ^ + E I i r Vl-I^MPe-^dw 

• /_7r rj£Z-{0} J ~ 7r 



Using 1 — Vl ^ ^ x, Vx e (0, 1), yl — x = 1 + YlpLi a P x p , that converges uniformly with 
respect to x £ [0,6] when <5 < 1 and where a p < for every p6N* and |39|) . we get 

N(z-i) < HZH^^j - £ «p E ^/:j^HI 2p e — ^ 

p=l n£Z-{0} 

oo 

*S iV(Z) 2 - £ a p 7V(|Z| 2 f) 
p=i 

OO 

N(Zf - £ a p N(Z) 2 P 
p=i 



iV(Z) 2 + 1 - y/1 -N(Z) 2 
< 2iV(Z) 2 . 

Formula (|4*2|) is a direct consequence of the previous inequalities. □ 
The goal of this section is the proof of the following theorem. 
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Theorem 7 There exists 8 > such that, for every Mq : [0, w] — > § 2 wii/i -/V[Zo] < S and 
zq < —1/2, there exists e = e(Mo) > swc/i i/iai, £/ie solution of (E) with M(0) = Mq, 



2fc-l 



(*) : — _ l[k,k+e](*) _ J! S ( C -fe+p(^o)) -l[fe+p,fe+p+e]( i ) + -l[3fc,3fc+e](*) 



p=l 
2fe-l 



P =i 

where k = fc(Mo) € N is such that 



N(Z ) 

|ra|>& 

satisfies 



E M*b)l < ^f 2 , (43) 

\n\>k 



N[Z(3k + e)]< N[Z ° ] 



2 ' 

z(3fc + e) < -1/2. 

By iterating this process, we find an increasing sequence (t n ) n eN € [0, +oo) N and two 
controls h,ij£ Lj£ c ([0, +oo),R) such that 

N[Z(t n )] < ^N[Z }. 

Thus, ||M(t n ) + C3||l°° — * when n — > +oo. These explicit controls provide the exact 
asymptotic controllability to e^. 

In section 15. 1[ we present the heuristic of the proof of Theorem [7J which is detailed in 
section 15.21 



5.1 Heuristic 

Let us sketch the proof of Theorem [7J It is inspired by the return method, introduced in 
[6] and already used for the control of quantum systems in 0] (for other applications 
see the book [7j). It consists here in going close to in order to delete the main Fourier 
coefficients of the initial condition, and then to move back to — e$. 
Notice that, when z > 0, the system ([3]) implies that 



Z(t,u) = iwZ(t,u) - w(t)y/l - \Z(t,u)\ 2 , (t, w) € (0,+oc) x (0,tt), (44) 

z(t,tj) = -U[w(t)Z(t,u>)],(t,u) e (0,+oo) x (0,tt). (45) 

We have 

Z (w) = E d « em "' where d » := c «( Z o)- (46) 

raGZ 

Let k GW that will be chosen later on. On the time interval [0, k) we take w = 0, thus 
Z(fc-,w) = Z (w)e* w = E^e l(n+fc) " and z(AT,w) = z (uj). 

At time k, we apply the control w(t) — inSk(t) in order to move close to +e%- Indeed, 
thanks to Definition [IJ we have 

1 

M(k + ,uj) = exp(Trn x )M(k-,uj) = \ 0-10 ) M{k~ ,w), 

0-1 
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thus 

Z{k+, u>) = Z(k~,uj) = d^e l{ - n ~ k)ul and z(k + , lu) = -z(k~,w). 
On the time interval (k, 3k) we apply a control of the form 

2&-1 

w(t) = ^2 Wp5 p+k (t), 
P =i 

where w p G C. Approaching the nonlinear system (|44ll by its linearized system around 
[Z = 0,w = 0), we get 



Z(3k~,w) ~ (z{k+ ,u) - !l k w^e-^V-VdtX 



_ 2k-i \ (47) 

,n(EZ p=l / 



Moreover, z stays close to +1 because the control applied is small. Choosing w p := d p -k, 
we get 

Z(3k-,uj) ~ ^ d^e l( - n+k ^. 

|n|>fe 

Finally, at time 3fc, we apply the control u>(i) = i7r<53fc(i) in order to return to —63: 
Z(3fc+,u) = Z(3fc-,w)~ 51 d «e l( "- fe) " 

and z(3/c + ,o;) = — z(3fc~,w) is close to —1. Now, by choosing k — k(Zo) such that 

\dn\ < \n{Z ), 

|n|>fc 

we get the existence of a time T — T{Zq) := 3k and a control w : [0, T] — > C such that 
JV[Z(T)] < AT[Z ]/2. 

Finally, the steps that need to be justified are 

• the approximation of the nonlinear system by its linearized system in (|47|) , 



the convergence, for the norm N, of the solutions of ((3]) when we approximate the 
Dirac controls by controls in Lf£ c . 



5.2 Proof of Theorem \7\ 

Let us recall that the solutions of |(3]) with Dirac controls have been defined in Definition [TJ 
and that we have the following result. 

Proposition 6 Let /3,7 G R, M G C°([0, tt], § 2 ) be such that N{M ) < +00. Let M be the 
solution of (dp with M(0) = M , u{t) = 0S o (t), and v(t) = jS (t). For e > 0, let M e be the 
(classical) solution of with M(0) = Mo, u(t) — (/3/e)l(o i£ )(t), and v(t) = (7/e)l(o,e)(i)- 
Then 

N(M e (e) - M(0+)) -> w/ien e -> 0. (48) 
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Proof of Proposition [6} We have 
M e (e,w) = exp[e\cj\fl z + f3fl x +jfl y }M (Lj), M(Q+,io) = exp[/3ft x +jQ y ]M (Lj). 
One has 

+00 , , 

N(M e (e) - M(0+)) < V (49) 
* — ' n! 

n—l 

with 

On(e) := N([(e\Lu\fl z + 0Q X + juj y ) n - (/30, + 7 ^)"]Af (^)) . (50) 

Noticing that 7V(|w|) < +00, using f39|) together with standard estimates and the Weier- 
strass M-test, one easily sees that (|4"8f follows from (|4"9|) and (f5Q|l . □ 

Thanks to Proposition [HI Theorem [7] is a consequence of the following result. 

Theorem 8 There exists S > smc/i iftai, /or every Mq : [0, 7r] — > § 2 with N[Zq] < 5 and 
z < -1/2, the solution of ® u»tft M(0) = M , 

2k-l 

u(t) := irS k (t) - ^ ^(c-k+p(Zo)jSk+p(t) + 7r<5 3fc (i), 
P =i 

2fe-l 



p=l 



where k = k(Mo) £ N is such that J^ffp holds, satisfies 



N[Z(3k+)} < ijV(Zo), (51) 

z(3fc+) < -1 (52) 
The key point of the proof of Theorem [8] is the following result. 

Proposition 7 There exist C > and C > such that, for every do <E C wii/i |<2q| ^ 1, for 
every Mq = (xo, 2/0, z o) : [0, 7r] — > S 2 wii/i iV(Zo) S% 1 and zo > 0, the solution of ^ w ^ 
M(0) = M , u(t) = -SR(d )*oW, "(0 = 3f(<2 )£o(*) safe/zes 

N(z{0+)-Z a + do) ^C\d \max{\d \,N(Z )}, (53) 

*(0+w)^«b(a;)-C / |do|max{|do| > JV(2b)}. (54) 
Proof of Proposition [fj Let us write do = Po + «7o ; with /?Oj7o G We have 
M(0 + ,w) = exp^o^x + 7o^]M (o;). 

Using the decomposition 

exp[/3 O x + 70^] =7 + /3 O a: +7ofi !/ +i2, where ||i?|| = O(\d \ 2 ) as d -► 0, (55) 
we get 

Z(0 + ,w) = Zo(w) - cfozoM + Rix (u}) + R 2 yo{uj) + i?3^o(w), 



24 



where Rj 6 C, \Rj\ < C | g?o 1 2 for j = 1,2,3, and C is a universal constant. Therefore, we 
have 



c 



„[Z(0 + ) - Z + d ] = d c n [l - zq] + RxCnlxo] + i? 2 c n [yo] + i?3C„[^o], Vn e Z. 



Using igOJl and we get 



JV(Z(0+) - Z + rfoj < |do|iV(z - 1) + |i?i|7V(x ) + |i? 2 |iV(y ) + |Jfe|JV(af ) 
< 2\d \N(Z ) 2 + C\d \ 2 [2N(Z ) + 1 + 2N(Z ) 2 }, 

which gives (J53J with C = 2 + 5C. 
From lf55"|) we get 

z(0 + , w) O^o(^)^ + R[x q (uj) + R' 2 y (uj) + R' 3 Zo(uj), 

where Rj 6 C, \R'j\ ^ C"|(io| 2 for j = 1,2,3, where C is another universal constant. Using 
(jlDjl and l[iT|l. we get 

«(0+,w) ^ z„M - MoPoHl ~ C'\d \ 2 [\x (uj)\ + \yo(u,)\ + \z (cu)\] 

> z (cj) - \do\N(Z ) - C'\do\ 2 [N(x ) + N(y ) + N{z )} 

> z (lo) - \d Q \N(Z Q ) - C'|do| 2 [2^(^o) + 1 + 2N(Z ) 2 ] 

which gives {P} with C = 1 + 5C".D 

Proof of Theorem [8t Let (5 be such that 

4C5<1, C'S< 1/2, 5 6 (0,1), (56) 

where C,C are as in Proposition [71 Let M , fc, u, w be as in Theorem[8l We use the notation 
(Hi. 

.Rrs£ step: on We have (see the previous section) 

Z(k + ,u) = 'fce ii ~ n ~ k)u and z(k + ,w) = -z (uj). (57) 



Second step: on (k,3k). Let us prove by induction on p £ {0, ...,2k — 1} that for every 
p 6 {0, 2fc — 1}, we have 



iv[z((fc+p)+)-E„ eZ - { _ fc+1 ,...,- fc+p} dne i (-"- fc+ ^l"l 

^c[|d_ fc+1 | + ... + |d_ fc+P |]^(z ), 



(58) 



(H' p ): z{{k+p) + ,u)>-z {u)-C'[\d^ k+1 \+... + \d^ k+p \}N{Z ). (59) 
Notice that {H^k-i) and ([56]) provide 



iV 



Z(3k~)- J2 d n e l{ - n+k) ^ ^CN(Z Q ) 2 , 



thus, thanks to l|56p and (|43"|) . we have 

JV[Z(3AT)] < JV[Z ]/2. (60) 
We also have, thanks to (-ff^fe-i) an d |56|) 

z(3k~,w) = z((3k-l) + ,u) > -z (w) +C'<5 2 > 
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thus 



z(3k~,u) = v/l - \Z(3k+,u)\ 2 > y/l-5/2 > 1/2. 



(61) 



The properties (H ) and (Hq) come from ([5?j) . Now, let p e {1, ...,2fc — 1} and let us 
assume that (H p _i) and (H' -y) hold. Thanks to (iJ p _i) and l|56p . we have 

+!>)")] = N[Z((k+p-l) + )] < JV[Z ] ^ <5< 1 
and thanks to (if^^) we have 

«((*+p)-,w) = z((fc+p- i)+ w) ^ -zqH - c'7V(z ) 2 > i-i = 

thus we can apply Proposition [3 Thanks to Proposition [7] and (iJ p _i), we get 
AT — • 



^ N 



Z({k+p)+) - d n e i( -- n - k+ P^ 
.ez-{-fe+i,...,-fe+p} 



?((fc + p)+) - Z((k + p)") + d_ fe+P 



Z((fc+p)- - £ d n e l (-"- fc+ ^ 

< C|<f_*+ p |JV[Z((* +?)")] 

+N\z((k+p- 1)+) - £ d^e^-™-^- 1 )^ 

- ^ ' nGZ-{-fc+l,...,-fc+p-l} 

s? C|cL fc+p |iV[Z ] + C[|d_ fc+ i| + ... + |cL fe+p _i|]7V(Z ), 



which proves (i? p ). Thanks to Proposition [7] and (H'l), we get 

z((k+p)+,u) > z({k + p)-,u;)-C , \d^ k+p \N[Z((k+p)-)] 
> z((k+p-l)+,uj)-C'\d- k+p \N[Z } 
^ -z Q {u) - C'[\d-k+i + •■• + |d-fc+p-i| + \d-k+ P \]N[Z ]. 



Third step: at 3k. We have Z(3k + ,uj) — Z(3k ,u>) and z(3k + ,uj) = — z(3k ,u>), thus 
((60)) and jSU) give ^ and (SSJ>. □ 



6 Comparison 

In this section, we compare the control results and processes presented in sections [U and [5j 

First, let us compare the statements of Theorems[6] (or Corollary 1) and[7l On one hand, 
the statement of Theorem[6]is stronger than the one of Theorem[7]because it is global and it 
gives the approximate controllability of $2$ for the norms ||.||#s,Vs < 1 (whereas Theorem 
[7] only provides the approximate controllability for N). On the other hand, Theorem [7] is 
stronger than Theorem [6] because it needs less regular initial data. 

Now, let us compare the control processes detailed in the proof of Theorems [6] and [7j 
Given Mq € H 1 ((oj*,li*),S ), the proofs of Lemma[T]and Proposition [5] give an explicit way 
to find T > 0, u, v, S D such that 

\\U[T+;u,v,M }'\\ L 2 < UQ L 2. 

Iterating this process, we produce a sequence of reachable points (Af„)„ S N c 
i/ 1 ((w >t ,u;*),§ 2 ) such that (j| Af^|| i 2)„ e N decreases. We expect that ||M^|| £ 2 — > when 
n — > +oo, and once this norm is small enough, we apply a control given in Lemma Q] (2) to 
go closer to e$. However, the sequence {M' n ) n£ m may not converge to 0. Thus, the control 
process presented in section [4] is not completely satisfying from a practical point of view. 
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Moreover, even if the sequence (M^)„ e m converges to in L 2 ((uj*,lj*) 7 M. 3 ), the control- 
lability process may take a long time (in particular the controllability time is not a priori 
bounded by a quantity depending only on ||Mo||jji) and cost a lot (because at each step, 
one has to compute new controls u and v and because the commands proposed in the proof 
of Proposition \E\ involve many trips between —63 and +63). 

On the contrary, the controllability process presented in section [5] works within a time T 
which is explicit, with controls u, v that are also explicit in terms of the Fourier coefficients 
of M , and needs only two trips between ±63. Thus, the time and the cost are well known. 

Let us compare the time and the cost involved by the two controllability processes on a 
particular example. We take (w*,cj*) = (0,7r/2), and an initial data of the form 

ex e (u>) 

Mo(cj) = I ' 



V'l - e 2 x e {ujf 



where e > is small, 

N 



x e (uj) = ^a fc (e)cos((2fc - + cos((27V + e (0,tt/2), (62) 



and (afc(e))i<jfe<gjv £ R are such that 



fc=i 

N 



tt/2 



f uj K duj = 0,VK G {0,...,JV-1}. (63) 



JO y/1 - € 2 X c (lo) 2 

We will prove later the existence of such coefficients. We want to reach e^. 

Let us apply the strategy presented in section HJ to find explicit T > 0, u, v S D such 
that 

\\U[T+;u,v,M }'\\ L 2 < \\M \\ L2 . 
One needs a polynomial Q S M[X] such that 

it/2 

(zx' — z'x)Qduj < 0. 

Then, deg(Q) > N, because of (|63|> . Thanks to the proof of Lemma [TJ there exists r* = 
r*(Q,x e ) > and a > such that, for every r 6 (0,r*), there exist T > 0, u,v £ D such 
that 

||(7[T+; U ,«,Afo]'|| 2 L2 < IIA^Hia-ar. 

However t* cannot be quantified, thus, we do not know the size of the decrease. Moreover, 
as emphasized in the proof of Proposition [U the time of control T satisfies T ^ 2 N t x / n 
(time needed to generate I + tu n fl x + o(t)) and one makes more than 2 N trips between ±e3 
(just count how many times the matrices exp(irfl x ) or exp(7rf2j,) appear in the generation 
of I + tuj n H x + o(t) in the proof of Proposition [5]) . 

With the strategy of section [5] taking the same explicit expression for Mq on (— w, n), we 
know the existence of e* > such that, for every e € (0, e*), the explicit controls 

u(t) := n52N+\{t) + n6 6N+3 (t), 

N 2N 

v(t) ■= — — ^ aN + i- m (e)82N+l+2m(t) — -z ^ am-Jv(e)^2JV+l+2m(*)! 
m — 1 m— N+l 
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with the convention a/v+i(e) = a-jv-i(e) = 1, realize 



TV 



U[{6N + 3)+;u,v,M } + e 3 ] < ~jV[M + e 3 ]. 



Here, the controls are explicit, the time scales like 6TV, we have a bound from below for the 
decrease of the TV-distance to and the process needs only 2 trips between ±e3. 

Now, let us prove the existence of the coefficients (cik(e))i^k^N ■ 

Lemma 3 Let TV G N* . 

(i) The matrix A G Mn(M) with coefficients 

,7T/2 

A k ,K ■■= / (2k - 1) sin((2A: - l)uj)uj K duj, l<fc<jV,0<.K'<jV-l, 
is invertible. 

(ii) There exists e* > and a C 1 map e £ [0, e*] i— > (afc(e))i^fc^w G such that 
EE-EM hold. 



Proof: (i) We assume that A is not invertible. Then, there exists (Ai, \n) SM^ — {0} 
such that 

/•tt/2 N 

/ VA fc sin((2/c- l)w)w K du = 0,V0 < K ^ TV- 1. (64) 

fc=i 

Let /(cj) := Ylk=i Afesin((2/c — l)w) and < w\ < ... < ojl < n/2 be all the values of 
the open interval (0,7r/2) on which / vanishes and changes its sign. Then, the function 
uj i — > f(u>)(u> — u>i)...(u> — lul) has a constant sign on (0,7r/2) and it is not identically zero, 
thus 

fTT/2 

f(u))(u> — U>i)...(u) — UJL,)dLJ ^ 0. 

The assumption (|64| ensures that L ^ TV. Thanks to trigonometric formulas, there exists 
(Hu • e E w - {0} such that 

N N 

/H = ^ ^fe sin(o;) 2fc - 1 = sin(c) £ ^ sinH 2 ^ 1 ) . 
fe=i fc=i 

Since the quantities sin^i) 2 , sin(u;jv) 2 are all different from zero (ui,...Uff 6 (0,7r/2)), 
they provide TV roots for the polynomial 

N 

1) 



E^ (fc_1) 



fc=l 



that have a degree ^ (TV — 1) and is different from zero. This is a contradiction, 
(ii) Thanks to (i), there exists (ai, ...,a/v) € such that 

«r/2 /n \ 
J ( X, ak ( 2k ~ ^ sin (( 2k - 1)^) + ( 2Ar + 1) sin((2TV + l)w)J lu k du = 0, V0 K < TV. 



There exists M > such that 

AT 



^a fc cos((2/c- l)w) + cos((2TV + l)w) s; Af,Vw G (0,tt/2). 



fe=i 
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When b = (b 1 ,...,b N ) t e R- 



we have 



N 



^(a fe + b k ) cos((2£; - l)w) + cos((2iV + l)w) < Af + ViV||6|| 



fe=i 



thus, the following map F is well defined 



F: (0,5*7) x S KN 




F(e,6) 




^1 - e 2 y b (w) 




where 



AT 



:= + &fc) cos((2fc - l)w) + cos((2iV + l)w), Vw e (0, tt/2). 



fe=i 

Then F(0,0) = and <4F(0,0) is invertible, thanks to (i). Thus, the implicit function 
theorem gives the conclusion. □ 
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